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Abstract 

We compute the cosmological background radiation of gamma rays and neutrinos 
due to neutralino annihilation in evolving dark matter halos, assuming the observed 
dark matter is comprised of thermally excited neutralinos in the MSSM. The spec- 
trum of this gamma-ray background radiation does not show strong annihilation 
line features, but could amount to a significant fraction of the extragalactic gamma 
ray continuum flux observed by EGRET above a few GeV. The corresponding cos- 
mological neutrino background is weak compared to the atmospheric foreground. 
Assuming full mixing, however, the cosmological tau-neutrino background could 
be detectable with a flavor-discriminating neutrino telescope in the energy range 10 
GeV - 1 TeV. A small anisotropy of the background radiation is expected, reflecting 
the local clustering of dark matter halos along the supergalactic plane. 
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1 Introduction 



The nature of the dark matter in the Universe is one of the biggest mysteries 
in modern astrophysics [1]. Of the numerous proposed constituents, the su- 
persymmetric WIMP (commonly believed to be the lightest neutralino Xi) is 
perhaps the most convincing candidate [2]. A possible way to find signatures 
of neutralino dark matter is to look for gamma rays being produced due to the 
self-annihilation of the neutralinos [3] . Detection of a gamma ray line signal at 
energies above a few GeV would constitute an immediate proof of new physics 
beyond the Standard Model. However, the process XX ~* 77 is loop- suppressed 
[4], and so these intensities are expected to be very small. Therefore, it has 
been proposed to search for continuum gamma rays resulting from the produc- 
tion and following decay of 7r-mesons (and to a lesser extent also K-mesons) in 
the annihilation process. Searches for continuum gamma rays originating from 
the dark halo of our own galaxy, but also from several extragalactic targets 
(including dSph-galaxies as well as M87, M31 and other nearby galaxies) have 
been discussed (see e.g. [5,6,7,8] for a review). Since neutralino annihilation 
has taken place throughout the Universe essentially since the freeze-out of the 
dark matter, a homogeneous and - at a first glance - isotropic background of 
annihilation products should also exist [9,10,11]. Comparing the measured ex- 
tragalactic gamma ray background and predictions from blazar models [12,13], 
one finds that there might be room for additional contributions. In this pa- 
per, we focus on the continuum gamma ray background (Section 2), as well 
as on expected high-energy neutrino backgrounds (Section 3). These annihi- 
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lation induced backgrounds are significant since they are the only cases of 
extragalactic backgrounds scaling with the matter density squared, instead 
of just the matter density, as for all conventional backgrounds. We discuss 
a possible anisotropy on large angular scales in those backgrounds owing to 
the anisotropic mass distribution in the local Universe, which is a potentially 
distinct signature of neutralino annihilation. Furthermore, in the case of su- 
persymmetric dark matter, these intensities are irreducible backgrounds for all 
observations, including targeted SUSY- searches with next-generation IACTs 
(Imaging Air shower Cerenkov Telescopes), which are operated in ON - OFF 
detection mode for background suppression. 

2 Gamma rays due to neutralino annihilations 

The differential gamma ray intensity per solid angle due to neutralino annihi- 
lations in any given dark matter distribution is 

* 7 (e) = ^ >< x J p 2 x XK[e(l+z),z}xdi [e(l + z)] ^ dz , (1) 

where df[e(l + z)] is the differential energy distribution of the produced 
gamma photons per annihilation event, (crv) is the self-annihilation cross sec- 
tion averaged over the thermal distribution of the dark matter particles, and 
m x is the particle mass, k [e (1 + z) , z] parameterizes the gamma ray attenua- 
tion on cosmological scales, and / (? ds with ds = ^ dz is the integral along 
the line of sight over the dark matter density squared. The factor 1/2 is in- 
cluded to compensate for evaluating twice each initial state going into (crv) 
[11]. While i(crv)/m x 2 and the differential energy distribution are fixed by 
the choice of the SUSY-model, the integral over the dark matter density de- 
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pends on the astrophysical abundance and distribution of the dark matter, 
specifically the value of Qdm and the dominant type of dark matter halo 
and subhalo profile. e = E/(z + l)is the observed energy of a photon emitted 
at energy E. To calculate the total cosmological intensity from dark matter 
annihilations, we start with the finite number of neutralinos produced dur- 
ing thermal freeze-out, which are depleted thereafter following the Boltzmann 
equation dn x /dt = — (ay) n 2 (1 + z) 3 . This approach has also been employed 
in [10] and yields 



^Wuualino , x ' , \ (°" V ) ^ 2 DM P^rit % 

^ {£) ~ 4^ X ^ > 

1 + z) 3 x /c[e(l + z),z] x T(z) 



J h (z) 

v ' 

In existing literature, most attention has been devoted to the cosmological 
annihilation into gamma ray lines, while we will take a different approach 
and focus on the characteristics and potentially detectable features of the 
extragalactic continuum signal. Therefore, some adaptations of the ingredients 
in Eq. 2 have to be included. On the " astrophysical side" , T (z) denotes the 
intensity multiplier function obtained from simulations of structure formation 
[14,15]. h(z) is given through 

h(z) = ^2 M (l + z) 3 + f2 K (l+z) 2 + fiA , (3) 

and for the "concordance" cosmology we will assume in the following that 
Q x ^ Q DM ^ 0.23, Q M ^ 0.27, Q K = and Q A = 0.73 [16]. For the Hubble- 
Parameter H we will work with the current best fit value of 
H = 71 km s _1 Mpc -1 [16]. p crit is 3H 2 /87rG , the overall (critical) density 
of the Universe for f2 = 1. As has been pointed out e.g. in [17], the annihila- 
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tion induced intensity heavily depends on the amount of structure in the dark 
matter. Instead of a constant "dumpiness factor", in this paper we include 
physics-motivated calculations of the z-dependence of the intensity multipli- 
cation function. For redshifts up to z = 20 for T (z) we will use calculations 
of the " dimensionless flux multiplier" from [14], which trace the process of 
structure and substructure formation. The results are in part based on the 
calculations and computer code presented in [15]. At higher redshifts, knowl- 
edge about structure in the Universe is rather limited, and we conservatively 
assume a completely smooth and structureless dark matter distribution at 
z > 20. One of the remaining questions is how much of the already formed ha- 
los are destroyed due to processes like tidal stripping and merger events. If this 
is the case for a sufficient number of halos, the relative intensity contribution 
due to annihilation events might actually peak quite early and then decline 
again, as indicated in [14]. However, as the authors acknowledge therein, this 
effect might well be overestimated in [14], e. g. because in N-body simulations 
the central regions of merging halos tend to survive. A more sophisticated 
treatment of the details of structure formation history should definitely be re- 
warding and perhaps settle this question. It should be noted that if this peak 
is in fact less prominent than assumed here, the intensity from the relatively 
local Universe will be even more significant with respect to the overall annihi- 
lation induced intensity. A lower mass cutoff of 10 6 solar masses for the dark 
halos is assumed, which represents the mass scale down to which the paradigm 
of hierarchical structure formation presently seems to be well established. We 
study two representative cases of dark matter halo profiles, the Moore et al. 
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profile [18] 




1.5 



(4) 



^scalc / \ ^scale 



and the NFW-type halo profile [19] 




(5) 



In both cases, the Sheth and Tormen (ST) mass function presented in [20] is 
assumed. We notice that in the recent literature, while some cases of (elliptical) 
galaxies with very shallow and low- mass dark halos have been reported [21], 
there has been convincing evidence of the validity of r^-like dark matter halo 
profiles from x-ray observations of Abell galaxy clusters, with (5 = 1.35 ± 0.21 
in the case of A2589 [22,23]. At least in the case of this cluster, the observed 
profile is not quite as steep as the Moore et al. profile, but has some excess mass 
with respect to the pure NFW case. As has also been suggested in [14], the 
Moore et al. case can therefore be considered an optimistic but not extreme 
scenario. For none of our neutralino models this scenario is constrained by 
expected gamma ray intensities from substructure in the Milky Way halo, if 
the external radius of the Milky Way clumps is taken to be their virial radius 
[25,11]. Expected extragalactic radio intensities from annihilation products of 
neutralinos in our models also are not larger than the intensities due to normal 
and radio galaxies [24]. Compared to the enhancement factors employed in 
previous literature (e.g. [10]), the calculations used here result in a nearly an 
order of magnitude larger present-day intensity enhancement of 1.5 x 10 7 for 
the Moore et al. halo profile. 
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For the cosmological distances of interest here, high energy gamma rays cannot 
propagate freely through the Universe owing to pair production with infrared- 
to-ultraviolet background photons. The details of gamma ray attenuation are 
very important for observations at cosmological distances, and have been sub- 
ject to scientific debate in recent years (e.g. [26,27,28]). The optical depth 
due to pair creation in the metagalactic radiation field (MRF) for a source at 
redshift z q and at the observed energy e is calculated to be 

Zq 2 OO ^ 

t 77 (£, z q ) = c J j J j-t^" x n ( z , E) x <7 77 (e, E, /i, z') dE d/i dz', (6) 
ooo 

where dl/dz' is the cosmological line element, 9 is the angle between inter- 
acting photons, /i = cos (6), n(z,E) is the MRF photon number density and 
<r 77 (e, E, //, z) the pair production cross section. The Fazio-Stecker- Relation 
presented in [28] allows us to estimate the redshift z cuto fr([£ (1 + z)] , z) at which 
the optical depth for a gamma photon of a given energy will be of the order 
one, so that we get 

k [e (1 + z) , z] = exp 

for the absorption of the neutralino-induced gamma rays. We use the best- 
fit-model Fazio Stecker Relation from [28]. Since reemission of photons with 
initial energies of some ten GeV occurs below 100 MeV, reemission will be 
neglected for our analysis of the high-energy gamma ray part of the spectrum. 
The high energy gamma rays are potentially most interesting because in this 
regime, the SUSY-contribution might have the highest significance against the 
steep conventional backgrounds. Since the calculations from [28] depend on 
the properties of the metagalactic radiation field, and thus on the star forma- 
tion history, they cannot reliably be extended beyond a redshift of roughly 



ZcutoffKl +z),z] 



(7) 
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five. For redshifts larger than five and up to z = 20, which for the gamma ray 
background we take as the upper limit z max of the integration in Eq. (2), we 
use more robust estimates for the gamma absorption at high redshifts from 
[29]. Due to the combined effect of redshift and structure formation, our end 
result is not really sensitive to the details of gamma ray attenuation at high 
redshifts. Now, having discussed the astrophysical constituents in Eq. (2), we 
turn to the SUSY-parameter set. Even in the Minimal Supersymmetric Stan- 
dard Model (MSSM), there are at least seven free parameters relevant to dark 
matter studies. We use the DarkSusy numerical package [30] to make large 
scans over the MSSM parameter space. The models are required to produce a 
neutralino not excluded by current experimental results [31]. Furthermore, the 
considered models are required to produce 0.025 < £\h 2 < 0.176. The higher 
limit is a little larger than the present best-fit value for Q DM [16] to include 
possible residual uncertainties in the concordance model. The lower limit ac- 
counts for the fact that the dark matter in principle need not be constituted 
entirely of neutralinos, but they should make at least a significant contribu- 
tion. For models that will not produce Q x = Odm we rescale the resulting 
intensities accordingly. The differential energy distribution for each model can 
also be calculated using DarkSusy. Fig. 1 shows a scatter plot of more than 
3000 neutralino models that match the criteria mentioned above. For each 
model we plot the neutralino rest mass and the parameter "R", which we de- 
fine as the ratio of the continuum gamma ray intensity created by the given 
model assuming the Moore et al. halo profile, and the extragalactic EGRET 
intensity (9 x 10~ 10 cm" 2 s" 1 sr" 1 GeV" 1 ) [32], at an energy of 30 GeV: 



(^Neutralino 

R = |30GeV- (8) 



total, EGRET 
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Fig. 1. R as defined in Section 2 for 3 x 10 3 cosmologically valid models plotted 
over neutralino rest mass. The Moore et al. halo profile case is assumed. The circle 
denotes the highly productive model discussed in detail in Section 2. 

At this energy, the contributions from highly productive models in the Moore 
et al. case will reach the ten-percent level, indicating that there might well be 
a noticeable component due to WIMP annihilation in the measured EGRB. 

In Fig. 2, for one model with m x = 672 GeV and (crv) = 6.1 x 10 -25 cm 3 s _1 
we plot the cosmological gamma ray signal from neutralino annihilations as 
a function of observed energy. For the measured extragalactic EGRET back- 
ground, we use the new determination by Strong et al. [33]. The bulk of the 
EGRB is likely due to unresolved, faint gamma ray emitting AGN (blazars) 
[12,13]. However, as Fig. 2 indicates, a dominant cosmological component 
at high energies could also be produced by cuspy halo profiles, a high self- 
annihilation cross section and a neutralino mass in the range of (500-900) GeV. 
As can be concluded from Fig. 1, for many of the calculated parameter sets 
and especially in the NFW case of halo profile, the neutralino induced intensity 
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will however not be the dominating component in the extragalactic gamma 
ray background. Thus, for singling out this signal from the conventional back- 
ground, highly selective criteria will have to be employed. Such a criterion can 
be a significant branching ratio into the gamma ray line channels. Almost in 
all cases however this line intensity will even be much smaller than the contin- 
uum intensities presented above, and identifying an asymmetrically broadened 
and very weak line signal at an unknown energy in many cases might not be 
possible. Still, there can be additional characteristics to the neutralino anni- 
hilation induced signal. In Fig. 2 the upper solid / dashed lines show the total 
neutralino annihilation induced background in our two cases of halo profile, 




Observed Gamma Ray Energy / GeV 



Fig. 2. Continuum gamma ray intensities in the Moore et al. / NFW case over 
observed energy for the neutralino model described in Section 2; the upper solid / 
dashed lines correspond to the total SUSY-induced intensity in the respective case, 
and the lower solid / dashed lines represent the estimated contribution from the local 
Universe (up to z = 0.01). The extragalactic EGRET intensity from [33,32] (stars / 
dots) is also shown for comparison. 



10 



and the lower solid / dashed lines correspond to the estimated contribution 
from the "local" (up to z = 0.01) Universe. This is done to point out that in 
case of a noticeable anisotropic dark matter distribution in the local Universe, 
this anisotropy will lead to an anisotropic extragalactic continuum intensity. 
From [34] we extract the spatial characteristics of the local Universe, and the 
locally produced gamma ray intensity from the direction of the "supergalac- 
tic plane" is rescaled with the overdensity <5 4 o = 0.48 for a supergalactic scale 
radius of 40 Mpc. Mainly due to redshift, this "foreground" becomes more dis- 
tinguished at higher energies, and for the highly productive model presented 
in Fig. 2 in the Moore et al. case it will reach approximately percent level 
respective to the total EGRB at 30 GeV (assuming the measured total ex- 
tragalactic EGRET intensity of 9 x lO" 10 cm" 2 s" 1 sr" 1 GeV" 1 at 30 GeV). 
Since the bulk of the measured EGRB is presently assumed to be due to more 
distant blazar sources, the anisotropy tracing the supergalactic plane might be 
a distinct feature of the intensity component produced by dark matter annihi- 
lations. In the NFW case the absolute intensities would be more than an order 
of magnitude lower, but the anisotropic behavior of the annihilation intensity 
itself is not sensitive to the details of the assumed profile. Such an anisotropy 
could in principle be detectable for next-generation satellite experiments (e.g. 
GLAST). Owing to the strong p 2 -scaling of the annihilation intensity, one 
would actually expect to "resolve" a significant part of the mass concentra- 
tions in the supergalactic plane, so that the discussed anisotropy should in 
fact be even more prominent than our calculations indicate. To study the ef- 
fects of gamma ray absorption, in Fig. 3 we plot the calculated extragalactic 
intensities from three fiducial neutralino models with masses ranging from 100 
GeV to 5 TeV. For each model we plot the relatively nearby (up to z = 0.03) 
and the " cosmological" (z larger than 0.03) contributions in comparison. The 
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Fig. 3. Calculated extragalactic gamma ray intensities for three fiducial neutralino 
models with (crv) = 1 x 1CP 26 cm 3 s _1 

and masses of m x = 100 GeV, 500 GeV and 5 TeV. The Moore et al. case of 
halo profile is assumed. 



intensity becomes increasingly suppressed with rising particle mass because 
of declining number density. For heavier neutralinos, gamma ray absorption 
becomes the limiting factor for the observed intensity, and for gamma energies 
above the 1000 GeV range, the intensity from the relatively nearby Universe 
actually exceeds the cosmological contribution. For annihilating WIMPs with 
masses of 1 TeV and beyond, this effect would make the proposed anisotropy 
even more significant. At energies beyond the ITeV-scale, the conventional 
extragalactic background from far-away sources is also expected to vanish due 
to gamma ray attenuation. 
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Observed Neutrino Energy / GeV 

Fig. 4. Cosmological r-neutrino intensity over observed energy for a trial neutralino 
model and the Moore et al. profile; we plot the overall (upper solid line) and the 
estimated relatively local (up to z = 0.01; lower solid line) contribution in compar- 
ison. Also shown are the upper bound on the r-neutrino intensity from hadronic 
sources (dotted line) [42,43] and the estimated r-neutrino intensity due to cosmic 
ray interactions in external galaxies [44] (dashed line). 

3 Neutrinos from neutralino annihilations 

The second very interesting channel is annihilation into high-energy neutri- 
nos. The calculations presented for the gamma rays in Section 2 also hold 
true for the neutrino channels, taking into account some modifications: The 
absorption function k for neutrinos can be set equal to unity for the energy 
range we are interested in. Also, the upper limit z max for the integration can 
be relaxed, and for our calculations was set to z max = 10 4 , since for greater 
redshifts a (9(100 GeV) neutrino is redshifted into the overwhelmingly strong 
stellar background below 10 MeV. Branching ratios into neutrino producing 
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channels show a wide range of possible values. They can approach order unity 
for bb-dominated channels, and can be of the order 0.1 for tt and also for 
W + W~ processes. The differential neutrino spectra for these processes also 
show some degree of variation. Generally speaking, spectra will be "soft" (i.e. 
with a shallow maximum at energies below one tenth of the neutralino mass) 
for the bb-dominated channel, and harder, with a sudden drop just below the 
X?-mass, for the W + W~-process [35,36]. 

A robust approach to estimate the total neutrino yield from hadronization 
and the following decay of charged 7r-mesons via n — > /i + and the following 
/i — > e + 2 v is straightforward. For the total hadron spectrum from the an- 
nihilation, one can use an approximation to the Hill spectrum [37,38], already 
taking into account that almost all produced hadrons are isospinsymmetric 
pions: 



with r] n = E 7r /m x . For the decay chains mentioned above, each of these charged 
pions will produce three neutrinos / antineutrinos, and the neutrino mul- 
tiplicity can be found by integration. The branching ratios and differential 
neutrino spectra can be obtained using DarkSusy, SusSpect [39] and the 
Pythia/JetSet Lund Monte Carlo [40] codes. Due to neutrino oscillations, 
the v e : : v T ratio of the neutrinos arriving on earth will be 
1:1:1 (max. mixing). In Fig. 4, we plot the r-neutrino background for a 
trial XX ~^ bb [35] - dominated model with (or) = 8.8 x 10~ 27 cm 3 s _1 and a 
neutralino rest mass of 230 GeV. The Moore et al. case of halo profile is as- 
sumed. Rescaling of the local contribution to take the supergalactic plane into 
account is done according to the procedure for the gamma rays. Again one can 



dN. 



7T 
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(9) 



dr], 
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see that the foreground contributes a small but potentially interesting part to 
the overall intensity. Since for electron and muon neutrinos, the atmospheric 
neutrino background will be approximately six orders of magnitude larger 
than our predicted intensity we calculate the r-neutrino intensity instead, for 
which the expected relevant backgrounds arise only from neutrino production 
in AGN sources and from cosmic ray interactions in external galaxies. It should 
be noted that the intensity from hadronic sources might well be weaker than 
the presented upper bound, making the annihilation induced intensity more 
significant. The cosmo logical neutrino background due to neutralino annihi- 
lations is therefore a new and interesting signature of the cold dark matter. 
In contrast to the well explored neutrino intensities arising from neutralino 
annihilations in the center of the Earth or the Sun, it is independent from the 
neutralino-nucleon scattering cross-section. If the neutralino is the WIMP, de- 
pending on the SUSY setup and strength of the AGN induced background, 
one would predict a characteristic "bump" in the extragalactic high-energy 
neutrino background at an energy corresponding to approximately ten per- 
cent of the neutralino mass. In future neutrino experiments, discrimination 
of r-neutrinos could be achieved by searching for a "double bang" signature 
(see e.g. [41]) in the detector, albeit rather in the time structure of the PMT 
response rather than in the spatial distribution of light along the track. 

4 Discussion 

From the results presented above several conclusions can be drawn: 

(1) For the neutralino models we calculated, the gamma ray background due 
to WIMP annihilation can constitute a sizable fraction of the total ex- 
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tragalactic gamma ray intensity. The anisotropies in the WIMP-induced 
background are expected to be different from those in the astrophysical 
background (from faint and distant gamma ray emitting sources) due to 
the unique p 2 -dependence of the dark matter induced intensity. 

(2) Owing to the expansion of the Universe and structure evolution, for both 
gamma rays and neutrinos a significant contribution to the overall an- 
nihilation induced intensity is produced in the relatively local Universe. 
The 7-ray absorption at high redshifts sharpens this effect, in particular 
for neutralinos with TeV-scale masses. 

(3) Taking into account that the mass distribution in the local Universe up to 
z 0.01 is not isotropic but concentrated in the so-called supergalactic 
plane [34], there should be a large angular-scale intensity anisotropy on 
the level of some percent (depending on the SUSY scenario and halo 
profiles). Small scale anisotropies due to nearby galaxies and dark matter 
clumps in the Milky Way halo should also exist [5,25]. 

(4) Using the halo profile and formalism given in [5], one can compare the 
estimated gamma ray intensity produced in the dark matter halo of the 
cD-galaxy M87 (a possible target for indirect SUSY-searches) and the cos- 
mological SUSY-induced intensity. We find that for a threshold energy of 
50 GeV within a radius of 6 arc minutes from the center of M87 the cos- 
mological intensity is about two orders of magnitude below the expected 
intensity from this individual halo. For targeted IACT observations, M87 
should stand out very clearly from this irreducible background. 

(5) High energy neutrino intensities from cosmological neutralino annihi- 
lations can be non-negligible when compared to other expected extra- 
galactic neutrino intensities [44], and should also show the discussed 
anisotropy. To make a detection of these WIMP-induced neutrinos feasi- 
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ble, one would need a flavor- discriminating z/ r -detector (see e.g. [45,46]) 
to suppress the very strong atmospheric u c + //^-background in the 10 
GeV - 1 TeV energy range. 

To summarize, neutralino self-annihilation is a potentially interesting pro- 
cess to produce cosmological gamma ray and neutrino backgrounds. As a 
corollary, these irreducible backgrounds must not be neglected in indirect 
WIMP searches. Further studies are demanding, in particular to incorporate 
mSUGRA-inspired regions of the parameter space into the scans and to rig- 
orously asses the density perturbation evolution. 
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